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I. INTRODUCTION 



The production of heavy quarks is important for our understanding of nature. Intensive 
experimental studies of a heavy quark production in various reactions involving unpolarized 
initial particles are presently being carried out. However, not less important are analogous 
studies when initial particles are polarized, either longitudinally or transversely. Experi- 
ments on longitudinally polarized initial particles are taking place [[), are being planned 
and proposed 0. These will deal with many aspects of polarized reactions, e.g. the prob- 
lem of the polarized gluon structure function g>2, the validity of the Drell-Hearn-Gerasimov 
(DHG) sum rule in QCD, etc... 

At the leading order (LO) Born term level, various heavy quark production mechanisms 
have been studied some time ago. However, the importance of knowing the next-to-leading 
order (NLO) corrections cannot be overemphasized. Unpolarized NLO corrections for heavy 
quark hadroproduction were first presented in , and in ||[7| for photoproduction. Cor- 
responding polarized results were calculated in || and |9,1C,11,T2|]. In all these papers 
cross sections were obtained by folding in the Born term matrix elements from the very 
beginning. Analytical results for the so called "virtual plus soft" terms were presented in 
for the photoproduction and unpolarized hadroproduction of heavy quarks. Com- 
plete analytic results for the polarized and unpolarized photoproduction, including real 
bremsstrahlung, can be found in [p^] . 

Let us emphasize the importance of knowing one-loop matrix elements, which contain 
the full spin information of the relevant subprocess. When the one-loop contributions are 
folded with the Born term contributions and spin summed, as in a NLO rate calculation, 
the information on the spin content of the one-loop contribution is lost and cannot be 
reconstructed from the rate expressions. On the other hand, having expressions for matrix 
elements allows one to easily derive the one-loop contributions to partonic cross section in- 
cluding any polarization of the incoming or outgoing particles. Also, it allows one to obtain 
any of the crossed processes, including the ones with a heavy incoming particle, that are 
needed for the different versions of variable flavor number schemes, for the direct checking 
of the DHG sum rule on the partonic level and possibly for doing future parametrizations 
for structure functions. And all this without doing calculations from the beginning. This 
work presents detailed results on a NLO calculation of partonic matrix elements for the 
set of one-loop Feynman graphs present in hadroproduction of heavy flavors, separately 
for every Feynman diagram in order to facilitate the use of the results for other relevant 
processes that differ by color factors. 

The one-loop results presented in this paper are an important input for part of the next- 
to-next-to-leading (NNLO) order calculation of the perturbative corrections to heavy flavor 
production. Last but not least, the imaginary parts of the one-loop amplitudes presented 
here are a necessary ingredient for the calculation of T-odd observables which are known 
to be fed by the NLO absorptive (or imaginary) parts of the one-loop contributions to a 
given process. 
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The subprocesses that receive one-loop corrections which are considered in this paper 
proceed through the following two partonic channels: 



9 + 9 -»• Q + Q, 



XV, 



where g denotes a gluon and Q(Q) denotes a heavy quark (antiquark), and 

q + q — > Q + Q, 



1.21 



where q(q) is a light massless quark (antiquark). 

We note that the Abelian part of the NLO result for Ql.lQ provides the NLO corrections 
to heavy flavor production by two on-shell photon collisions 



;i.s) 



with the appropriate color factor substitutions. The results for can be also used to 
determine corresponding amplitudes for heavy flavor photoproduction 



i + g -> Q + Q- 



;i.4) 



We mention that the partonic processes ( |1.1| ) and (O) are needed for the calculation of the 
contributions of single- and double-resolved photons in the photonic processes (|1.3|) and 
(0)- 

Cross sections for the process 



]1.3|) have already been determined in p3| ; [14| , [15| for un- 
polarized and in |T^|,|TBJ for polarized initial photons. Note that the authors of HJ used a 
nondimensional regularization scheme to regularize the poles of divergent integrals. In the 
papers ||T3|,|TB| analytic results were presented for "virtual plus soft" contributions alone. 
We also note that complete analytic results including hard gluon contributions can be 
found only in |3j|]. The reaction (|1.3|) will be investigated at future linear colliders. NLO 



corrections for the heavy quark production cross section (|1.3| ) are of interest in themselves 
as they represent an irreducible background to the intermediate Higgs boson searches for 
Higgs masses in the range of 90 to 160 GeV. 

The paper is organized as follows. Section II contains an outline of our general ap- 
proach as well as matrix elements of the gluon fusion subprocess for the self-energy and 
vertex contributions including their renormalization. In Section III we discuss the one- 
loop contributions to the four box diagrams in the same gluon-gluon subprocess and give 
a detailed description of our global checks on gauge invariance for our results. Section IV 
presents analytic results on the quark-antiquark subprocess (|1.2j ). Whereas Sections II-IV 
deal with the real (or dispersive) parts of the one-loop amplitudes we turn to discuss their 
imaginary (or absorptive) parts in Section V. Section V also includes a discussion on how 
one may obtain the corresponding absorptive parts in different crossed channels. Our main 
results are summarized in Section VI. Finally, in two Appendices we present results for 
various coefficient functions completing our determination of NLO matrix elements. 
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II. CONTRIBUTIONS OF THE TWO- AND THREE-POINT 
FUNCTIONS TO GLUON FUSION 



The Born and the one-loop contributions to the gluon fusion partonic reaction g(pi) + 
g{P2) Qipz) + QiPi) are shown in Figs. 1-3. In this section we discuss our calculation 
of the self-energy and vertex graphs that contribute to the above subprocess. With the 
4-momenta pi{i = 1, ...,4) as indicated in the Fig. 1 and with m the heavy quark mass we 
define: 

s = (pi+p 2 ) 2 1 t = T -m 2 = (p 1 -p 3 ) 2 -m 2 , u = U-m 2 = (p 2 -pz) 2 -m 2 . (2.1) 

To isolate ultraviolet (UV) and infrared/collinear (IR/M) divergences we have carried 
out all our calculations in both conventional regularization schemes, namely the standard 
dimensional regularization scheme (DREG) |T7| and the dimensional reduction scheme 
(DRED) [Tj|. In what follows, we present results for the DREG, as well as the difference 
A=DRED-DREG. A brief characterization of the two regularization schemes is the follow- 
ing: In DREG both tensorial structures (e.g. gamma matrices, metric tensors, etc..) and 
momenta are continued to n ^ 4, while in DRED only momenta are continued to n ^ 4 
whereas the tensorial structures are those of n — 4. 

First of all we note that in general the matrix elements for all the Feynman diagrams 
in the gluon fusion subprocess are written in the form 



However, for the purposes of brevity, we will present our results in terms of the amplitudes 
M^ v omitting the polarization vectors and Dirac spinors. Of course, their presence is 
implicitly understood throughout this paper in that the mass shell conditions Pi^^pi) = 
and p^ 3 u{p 3 ) = mu{ps) etc. are being used to simplify M^ v . Furthermore, M^ u for all 
the one-loop graphs considered in this paper contains a common factor due to the one-loop 
integration. Define the quantity 



where g is the renormalized coupling constant. 

1 According to the discussion in ]i"9|] this implies that, when further processing our LO and one- loop 
results in cross section calculations by folding in the appropriate amplitudes, one may use the Feynman 
gauge for the spin sums of polarization vectors. At the same time ghost contributions associated with 
external gluons have to be omitted. 



(2.2) 




(2.3) 



We will omit from all of our one-loop amplitudes the common factor 



C = g i C £ {m 2 ) 



(2.4) 
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For our analysis of matrix elements it is important to describe various crossed heavy 
flavor production channels. We make it clear from the outset that an additional u-channel 
set of graphs, that topologically differ from the t-channel ones, are obtained by interchange 
of bosonic lines (not momenta). In particular, for calculational purposes, we will always be 
relating t- and u-channel Feynman diagrams by the following procedure: 

M t <-> M u = {a <-> b, pi^p 2 , li <-> v}, (2.5) 

with a, b color indices of bosons and where all three interchanges are performed simultane- 
ously. Note that the second interchange in Q2.5|) implies also the interchange t <-> u, but 
not vice versa. One case, involving two vertex diagrams, when the above transformation 
(|2.5|) does not correspond to "true" u-channel topologies, is discussed below. In general, 
when speaking about t-u symmetry of given amplitudes, we will imply invariance of those 
amplitudes under the transformations (|2.5| ). 

We start by writing down matrix elements for the leading order Born terms. For the 
t-channel gluon fusion subprocess (first graph in Fig. 1) we have: 

Br = -tT b T a r(4>3 "A + m)>f/t, 

where T b and T a are generators (T a = A a /2, a = 1, 8 and the A a are Gell-Mann matrices) 
that define the fundamental representation of the Lie algebra of the color SU(3) group. 
Analogously, for the u- and s-channels we have, respectively, 

= -tT a T b Y(J> 3 -fo + m)Y/u, 
Bf = i(T a T b -T b T a )Cr a la/s, 

where the tensor C% ua is defined according to the Feynman rules for the three-gluon cou- 
pling. We have omitted a common factor g 2 in the Born amplitudes. Acting with Dirac 
spinors on the above Born matrix elements from the left and the right and using the effec- 
tive relations p± = p v 2 = 0, as remarked on before, we arrive at the following expressions 
for the leading order matrix elements: 

BT = iT b T a { 1 ^ ll »-2pW)/t, 

Bj? = iT a T b (2pW-Yhin/u = iT a T b (Yfal ll -^)/u- (2.6) 
B? = 2i(T a T b -T^Wfa+pZY -plY)/s- 

Next we proceed with the description of the two-point contributions to the matrix 
element of the subprocess ( |1.1|) . But before we turn to the two-point functions one should 
mention that our choice of renormalization scheme will be a fixed flavor scheme throughout 
this paper. This implies that we have a total number of flavors rif = riif + 1, where is 
the number of light (e.g. massless) flavors plus one produced heavy flavor, with only nif 
light flavors involved/active in the /3 function for the running a QCD coupling a s and in the 
splitting functions that determine the evolution of the structure functions. When having 
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massless particles in the loops we are using the standard MS scheme, while the contribution 
of a heavy quark loop in the gluon self-energy with on-shell external legs is subtracted out 
entirely. 

Consider first the two t-channel self-energy graphs (2d2) and (2d3) with external legs 
on-shell (note that in the graph numeration the first number identifies the figure which 
the given diagram refers to). These graphs are very important as they determine the 
renormalization parameters in the quark sector. Throughout this paper we use the so called 
on-shell prescription for the renormalization of heavy quarks, the essential ingredients of 
which we describe in the following. When dealing with massive quarks one has to choose a 
parameter to which one renormalizes the heavy quark mass. It is natural to choose a quark 
pole mass for such a parameter - the only "stable" mass parameter in QCD. The condition 
on the renormalized heavy quark self-energy S r (j6) is 

£ r (i>)li/= m = 0, (2.7) 

which removes the singular internal propagator in these self-energy diagrams. The above 
condition determines the mass renormalization constant Z m . For the wave function renor- 
malization we have used the usual condition (see e.g. Ref. ||) 

^E r (^ =m = 0, (2.8) 

which fully determines the wave function renormalization constant Z 2 . Since the condition 
( |2.8| ) is not mandatory in general, there is a freedom in determining the constant Z 2 . 
Therefore, we will list our expressions for these constants. In DREG we arrive at 

Z m = 1 - 3g 2 C F C £ (m 2 ) (~ + ~) , Z 2 = l- g 2 C F C £ (m 2 ) (± + 4 + £) . (2.9) 



2 > / 1 _ 2 



And in DRED we obtain 

Z m =l- 3g 2 C F C £ {m 2 ) (i + ~) , Z 2 = 1 - g 2 C F C £ (m 2 ) (J + 5 + , (2.10) 

where C F =4/3 and where we use 1/e' to indicate which terms are of ultraviolet origin. Now 
we are in a position to write our results for the (2d2) and (2d3) quark self-energy diagrams. 
In DREG we have: 

Mfa) = Mfa = ~C F Br (p + 4 In , (2.11) 

where we clearly see the soft divergence in the logarithm that diverges on the mass-shell. 
Apparently, this very same logarithm appears also when calculating the wave function 
renormalization constant Z 2 . This means that on-shell mass renormalization does not 
actually completely remove the divergent propagator, but rather transforms the strong 
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linear divergence to the "softer" logarithmic divergence, which can be dealt with in QCD. 
Indeed, one can rigoriously prove that there exists the following effective correspondence 
relation: 

m 2 — p 2 1 1 1 

— — | p2=m2 ^ - + = - 

With this in mind we obtain final results for the two self-energy graphs after mass renor- 
malization was carried out: 



In | p2=m2 ->— + - - = — + i + 0(e). (2.12) 



DREG : M ( 7 d2) = M ( 7 d3) = -C F BjT ( 1 + 4 + ~ ) , (2.13) 



£• e 

DRED : M ( 7 d2) = Mf/ d3) = -C^f (p + 5 + |) , (2.14) 

and the difference between the two regularization schemes is 

A(2d2) = A(2d3) = -CpB^. (2.15) 

One notices that the effect of the wave function renormalization consists of a complete 
removal of the quark self-energy diagrams with external legs on-shell, as it is required by 
the second condition (|2.8|) . We can also write the contribution of the quark self-energy with 
external legs off-shell, graph (2dl), after addition of the mass renormalization counterterm: 

M (Tdi) = CfBF (-l/e' -t/T + \n(-t/m 2 )(4t/T + t 2 /T 2 -4)) (2.16) 
- iC F T b T a mYY (l - 2 ln(-t/m 2 ) - ln{-t/m 2 )t/T) /T. 

The difference between the DRED and DREG results is 

A(2dl) = -C F B?. (2.17) 

The remaining quark self-energy diagrams (3il) and (3i2) with external on-shell legs 
are derived in analogy to the ones considered above: 

Mfa = M$ 2) = -C F B? (1 + 4 + , (2.18) 

A(3il) = A(3i2) = -C F B?. (2.19) 

Concerning the gluon self-energy graphs (2el) and (2e2) with external legs on-shell, the 
only nonvanishing contribution they receive from the loop with internal heavy quarks is 
given by 

Mfci) = M&> = ~^ V \ \ (2-20) 

However, these contributions are explicitly subtracted (together with the logarithmic term 
ln(/i 2 /m 2 ) coming from the common factor C £ (m 2 ), see eqs. ( |2.3|) and Q2.4Q ) in the on-shell 
renormalization prescription, in order to avoid the appearance of the large mass logarithms 
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from the gluon self-energy diagrams with off-shell external legs in the low energy limit. 
Therefore, due to the UV counterterm that subtracts that very same loop with heavy 
quarks, there are no finite contributions to the matrix element from these diagrams. How- 
ever, at the same time this counterterm introduces the pole terms from the light quark 
loop sector that are needed to cancel soft and collinear poles from the other parts of the 
amplitude, e.g. from the real bremsstrahlung part. This indicates that in practice it is very 
hard to completely disentangle UV and IR/M poles in heavy flavor production and in most 
cases one obtains a mixture of both instead. 

For the reasons specified above it is convenient to present a gauge field renormalization 
constant Z 3 , used for the gluon self-energy subtraction: 



i + 4 



-ni f )C £ (^ 2 ) 



(b-2N c )C £ (/i 2 



- 2 c 

\c e {m 2 ) 



m 



(2.21) 



with the QCD beta function b = (lliVc — 2n^)/6 containing only light quarks. Nc = 3 is 
the number of colors. 

Similarly to the diagrams (2el) and (2e2), diagrams (3jl) and (3j2) also vanish due 
to the explicit decoupling of the heavy quarks in our subtraction prescription. However, 
instead of doing renormalizations separately for each Feynman diagram, one can chose to 
employ the renormalization group invariance of the cross section and do only a mass and 
coupling constant renormalization. In that case, knowing the results for gluon self-energies 
turns out to be useful in checking the complete cancellation of UV poles by just rescaling 
the coupling constant in the LO terms (?bare — > Z g g: 

-5f - -. (2.22) 



M (3j2) 



Finally we arrive at the gluon-self energy graph (3h), which contains the off-shell gluon 
self-energy loop that is used for the derivation of the renormalization constant Z%. We have 
evaluated the internal loop in the Feynman gauge. Since it is explicitly gauge invariant, we 
should arrive at the same result in any other gauge. In our result we show separately the 
gauge invariant pieces for gluon plus ghost, light quarks and a heavy quark flow inside the 
loop: 



M (3h) 



B 



A , ,15 31\ (1 2 10\ 

Nc[- o + — ) -ni f -7 



9 



9 J 



--i 



CJm 2 ) e' 3 



(2.23) 



with 



'5 Am 2 / 2m 2 \ fll . .1 
- + — + 1 + — W) . 



In fl2.24j ) we have made use of the definitions 

P = Jl-4m 2 /s, 



x 



1 + /?' 



(2.24) 



(2.25) 
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where (3 is the velocity of the heavy quark. 

We emphasize that in the last term of ( |2.23|) the Z 3 counterterm together with the UV 



pole will remove also the ln(/i 2 /m 2 ) contribution, while ln(— /i 2 /s) from the first two terms 
in (|2.23|) will be left unsubtracted. Expression ( |2.23| ) is obtained for the physical condition 



s > 4m 2 for producing a heavy quark pair. However, if one is interested in the case when 
s < Am 2 , then in fl2j| ) one should make the replacement 

(3ln(x) -2y / 4m 2 /s - 1 arctan ~ T = 1 (2.26) 

J Am? / s — 1 

We note that there is a minor problem with preserving gauge invariance when calculating 
the graph (3h) in DRED. It is associated with an e- dimensional part of one of the n- 
dimensional metric tensors g>" that arises in every partonic loop and hampers collecting 
together similar terms. However, this problem appears to be an artificial one, as in this 
particular case it makes no difference whether one uses 4- or n-dimensional metric tensor for 
the evaluation of this gluon self-energy graph. For this reason in practice one would set this 
g™^ metric tensor to be the 4-dimensional one. Or, more exactly, if one introduces a proper 
counterterm so that to restore gauge invariance of the gluon self-energy, then the expression 
in DRED would be exactly the same as in ( [2.23 ), except for the term proportional to Nc, 



where one would have 28/9 instead of 31/9, leading to 

A{3h) = -B%\ (2.27) 

Concluding our discussion on the 2-point functions we remark that the matrix elements 
for the additional u-channel 2-point functions can be obtained from eqs. ( |2.13|) , ( |2.16|) and 
(|2.20|) by the transformation (|2.5|) . 



We start by considering the t- and u-channel vertex diagrams. In particular we begin 
with the purely nonabelian graph (2b), which contains a four-point gluon vertex. It is 
finite, e.g. does not have UV and IR poles. For convenience we define the function 

/ Ud = 4Li 2 (-x) + ln 2 (:r) + 2((2), (2.28) 

with C(2) = *£. 

Then the matrix element takes the form 

^(Tb) = iN c {T b T a [A{2mg^ - 2p^ ~vVi» + PsY + 2pW){sf3\n{s/m 2 ) - m 2 f lld ) - 
3msYl u (4l3\n(s/m 2 ) - f m )}/As 2 (3 3 + (a <-> b, /i <-► v)} - 
i5 ab [2{2mg^ + p^ - p^ + p^ - pfr") {s(3 ln( S /m 2 ) - m 2 f m ) - 
3msf3 2 g^f m ]/As 2 (3 3 . (2.29) 

It is easily seen from eq. Q2.29D that graph (2b) is explicitly t-u symmetric, as it follows from 
its geometric topology. Because M/^\ does not exibit any UV poles, there is no difference 
between DREG and DRED in this case, e.g. 

A(26) = 0. (2.30) 
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Next we turn to graphs (2cl) and (2c2) and define another useful function: 

/ dl = C(2)-Li 2 (^). (2.31) 

Diagrams of this topology do not only occur in hadropro duct ion, but also in other processes 
such as photoproduction and 77 production of heavy flavors. For this reason we also present 
the corresponding t-channel color factors for these graphs. Then it is straightforward to 
separate our Dirac structure from the color coefficients and one can easily deduce results 
for other processes where these graphs contribute, though with different color weights. The 
color factor for both (2cl) and (2c2) diagrams is the same: 

T c (2cl) = T c (2c2) = (C F - ^)T b T a = --T b T a . (2.32) 

2 6 



The complete matrix elements are: 

fflV 

'(2d) 



MZi) = Bri-l/e' + 2f dl m 2 /t-\n(-t/m 2 )(Qm 2 +t)/T]/6 



+ iT b T a {p%Y [2f d m 2 T/t 2 + ln{-t/m 2 ){m 2 /T - 2T/t) + 1] (2.33) 
- mp%fii"/ v [\n(-t/m 2 )/T - 1/t) + m-f^ \n(-t/m 2 )}/3T. 

Because one has a UV pole 1/e' (the prime denotes UV poles), there is a finite difference 
between the two n-dimensional schemes: 

A(2cl) = (CV - ^)IT • (2.34) 

As expected, the difference is proportional to the sum of the Born terms of the relevant pro- 
duction channels. We emphasize that one has to take care when calculating matrix elements 
in DRED. In particular, one has to keep a clear distinction between the 4-dimensional and 
the n-dimensional metric tensors g^ v when doing their convolutions. Failing to do so may 
very well result in introducing extra terms that do not satisfy the Slavnov- Taylor identities 

For the graph (2c2) we obtain: 

^(Tc2) = Br [-!/£' + 2/ dl m 2 /t-ln(-t/m 2 )(6m 2 + t)/T]/6 

+ iT b T a {plY[-2fdirn 2 T/t 2 + \n(-t/m 2 )(m 2 /T + 2T/t)-2m 2 /t-l] (2.35) 

+ m^7 / ^ 1 [ln(-t/m 2 )/T - 1/t] + m^'-f ln(— t/m ) 

- 2mp%p u 4 [\n(-t/m 2 )/T - l/t]}/3T 



and 



A(2c2) = (C F - ?y)Br. (2.36) 
To write down the results for graphs (2c3) and (2c4) we introduce one more function: 

/n = ln 2 (-^) + Li 2 (^). (2.37) 
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The color factors for both diagrams are the same: 



rp(2c3) _ rp{2c4) _ jVg rpbrpa _ _^_J^brpa ^2 

Graphs (2c3) and (2c4) exhibit a rich structure of UV and IR singularities. Here one of 
the scalar integrals that is needed for the reduction of tensor integrals is of the type 



/ 



d k , n x 

(2.39) 



P(fc + p) 2 ' 



which is effectively zero in dimensional regularization. However, we choose to separate 
UV and IR/M poles to keep trace of all the sources of UV singularities, e.g. the above 
mentioned integral is proportional to the difference 1/e' — 1/e. This is our procedure for 
the relevant vertex diagrams. Thus, 

M (2c3) = 3Br[3/e' - 1/e 2 + (2 ln(-t/m 2 ) - l)/e + 4 + 6 \n(-t/m 2 )m 2 /T - 2/u]/2 

+ 3iT b T a { P %Y{-l/e 2 + 2/£ + 21n(-t/m 2 )/£-2/ u + 21n(-t/m 2 ) (2.40) 

x m 2 (m 2 /T+ l)/T + 2m 2 /T + 6]/2t - 3m^Y ln(-t/m 2 )/2T 

- mPs A7^[ln(-t/m 2 )(m 2 + T) + T]/tT 2 } 



with 



And 



A(2c3) = —Br. (2.41) 
2 



m qL) = 3ST[3/e / - 1/e 2 + (21n(-t/m 2 ) - l)/e + 4 + 6 ln(-t/m 2 )m 2 /T - 2/ u ]/2 

+ 3zT fe T a K 7 ^[l/£ 2 -2/e-2 ln{-t/m 2 )/e + 2/ n + 2 ln(-t/m 2 ) (2.42) 

x m 2 (m 2 /T+ l)/T + 2m 2 /T-6]/2t - 3m^Y ln(-t/m 2 )/2T 

+ m^7^ 1 [ln(-t/m 2 )(m 2 + T) + T]/tT 2 - 2mp^[ln(-t/m 2 )(m 2 + T) + T]/tT 2 } 

with 

A(2c4) = ^^r. (2.43) 

The results for the matrix elements of the additional u-channel vertex graphs are ob- 
tained from eqs. Q2.33| ),( f2T3l)| ),( |2.40"D and Q2.42p by the transformation (|2.5p. However, there 



is a subtle point involved here: we stress that for the graphs (2c3) and (2c4) transforma- 
tion Q2.5D transforms the t-channel result of the graph (2c3) to the u-channel result for the 
graph (2c4), while the t-channel result of (2c4) goes to the u-channel result for (2c3). This 
is important to keep in mind when dealing with reactions which involve asymmetric set of 
graphs as e.g. photoproduction of heavy flavors. 

It is worthwhile to mention that we have tested the above results for the t- and u- 
channel contributions against the ones given in |E|] and [ITH] by folding our NLO matrix 
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elements with the Born term matrix elements for the cases of unpolarized and longitudinally 
polarized incoming bosons. We obtain full agreement. 

Next we turn to the remaining s-channel graphs shown in Fig. 3. For all the gluon 
propagators we work in Feynman gauge. Although this set of graphs is purely nonabelian for 
QCD type one-loop corrections, there could be also abelian (e.g. QED) virtual corrections 
to graph (3fl). For this reason we also give the color factor for it separately: 

T c ( 3 ! fl) = (C F - ^-)(T a T b - T b T a ) = --(T a T b - T b T a ). (2.44) 

In the result for this graph together with the UV divergence we clearly see a collinear pole as 
well, multiplied by a logarithmic factor, that should be cancelled against the corresponding 
term in the factorization counterterm of the real bremsstrahlung: 

M^ {1) = {-B^s[s(3/e' + (2m 2 - s)(2\n(x)/e - ALi 2 (x) - 41m» ln(l - x) + 61n(x) 

+ ln 2 (a;) -8C(2)) + 31n(a;)s] (2.45) 
+ 2i(T a T b - T b T a )mln{x)[-g^{s + 2t) - Ap&Z + 4p^]}/6s 2 /3. 

A( 3 fi) = (C F - ^)B?. (2-46) 

Graph (3f2) contributes as: 

Km = N c {B^[3sf3 2 (l/e' + 2)-\n(s/m 2 )(8m 2 -s) + f m m 2 /f3] 

+ 2i(T a T b -T b T a )m[g» u (s + 2t)+4p%p u 4 -4p%p u 3 } (2.47) 
x [ln(s/m 2 )(8m 2 + s)/f3 2 - 2s - 3f m m 2 /(3 3 }/s 2 }/2s(3 2 . 

The difference between the two regularization schemes is again determined by the coefficient 
of the UV pole: 

A (3 f 2 ) = • (2-48) 

We finish our consideration of the vertex diagrams for gluon fusion with the triangle 
graph contribution (tri) = (3gl) + (3g2), e.g. we sum the two graphs (3gl) and (3g2). For 
the case when one has gluons and ghosts inside the triangle loop we obtain: 

M^ig) = N c {Bf[ZZ/ef - 36/e 2 - 171/e + 361n(s/m 2 )/e + 1381n(s/m 2 ) - 181n 2 (s/m 2 ) 
+ 144C(2) - 284] + Qi(T a T b - T^)^ [g^{21/e' - Q/e 2 - 33/ e + 6 \n(s/m 2 )/e 
+ 61n(s/m 2 ) -31n 2 (s/m 2 ) + 24((2) - 4)/s - p%p u 1 lQ/s 2 }}/72; (2.49) 

For the two more cases when one has light and heavy quarks inside the loop we get 

Mfab) = 2n lf {B^[3/e' - 31n( S /m 2 ) + 5] - 3t(T a T b - T b T a )^ '/ 's - 2p^/s 2 ]}/9 

(2.50) 
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with riif number of light flavors in the triangle loop, while for the heavy flavor case one has 

M^iQ) = 2{5f [3/e' + 31n(x)(2m 2 /s+l)/? + 5 + 12m 2 /s] (2.51) 

- 3i(T a T b - T h T a )i> x \g^ '/ 1 s 2 - 2p^/s 3 ] [3(ln(s) + 4/3) ln(x)m 2 

- 18C(2)m 2 + 24m 2 + s] } /9. 

The complete matrix element for the triangle is the sum of the above three expressions 
(g|9[) , (p30|) and flTSU ): 

^S) = M Sn)(^) + M &(<?) + (2-52) 

We have compared our results in eqs. ( |2.49|) and ( [2.50D with those available in the literature 
| f?T|1 , and found agreement. 

The difference between dimensional reduction and regularization arises due to gluons 
in the loop: 

A (tri) = Bf. (2.53) 



III. RESULTS FOR THE BOX DIAGRAMS IN GLUON 

FUSION 

In this chapter we describe the technically most complicated derivation of the 4-point 
massive loop diagrams. The results of this chapter cannot be deduced from the results of 
any other relevant publications up to date. The four box graphs (2al)-(2a4) contributing 
to the subprocess g+g — ► Q + Q are depicted in Fig. 2. We have used an adapted version of 
the Passarino-Veltman techniques [^] to reduce tensor integrals to scalar ones. The scalar 
integrals are taken from || which have been checked by us and in ||. 

First of all we note that the results for the box diagrams are the same both in DREG 
and DRED, as a consequence of the ultraviolet convergence of the box graphs. 

We expand all the box diagrams arising at the one loop level in gluon fusion in terms 
of eight independent Dirac structures. In turn, coefficients of the Dirac structures are 
expanded in terms of universal independent Lorentz objects. And finally, the coefficients 
of the Lorentz structures are expanded as products of a small set of analytic functions and 
various coefficient functions, which are combinations of scalar Mandelstam variables of the 
subprocess under consideration. Note that the Dirac structures and Lorentz objects are 
the same for every box diagram. In particular, we cast the box matrix element into the 
following universal form: 

= zT col {M^E/A (3.1) 

+ A br E ft® + m E /A? + pgrf E ft® + v>i E /4? + m E fJ$) 

+ 7>3 E + Pi E /#] + 7>£ E ft® +P4E /#] + mYY E fi<H 
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+ m^JPi [P3 fidn M2] + m 7% [P3 I] Mi + P4 H fea\ 

+ I] + P3P3 H /^*2 + P3P4 fi9is + P4P3 £ /i<&4 + £ /ife] } 

Even though the number of independent covariants for the process g + g—>-Q + Q is 
eight for n = 4 this will not be the case for 11 ^ 4 relevant to this application. We have 
therefore made no attempt to reduce the above number of covariants to a minimal set. 

Depending on the type of the box graph one has different number of terms under the 



summation signs in (|3.1|) . These numbers as well as the analytic functions fi are specified 
below. The coefficient functions bi, ...,^5 are given in Appendix A of this paper. 

The t-channel Born term matrix structure is the same for all box graphs and is 
defined as 

M^=Y(p 3 -Pi + m)Y, (3.2) 

which, when taken between the spin wave functions implying the effective relations p± = 
0, P2 = 0, can be written as 

M£ = 2p£7" - YhY- (3-3) 

Furthermore, for the box diagrams (2al) and (2a2) we empirically found the following 
relations between the d and e coefficients: 

e,;i = d i2 , e i2 = dn. (3.4) 

For the 4-point graph (2a4) the above relations are slightly changed: 

en = d i2 , e i2 = dii, for £ = 1,2,4,5,8,9; 

e3i = d 32 + 2m 2 u/(m 2 s -ut) 2 , e 6i = 2d 62 , e n = 2d 72 , (3.5) 

e32 = d 3l -2m 2 t/(m 2 s - ut) 2 , e 62 = 2d 61 , e 72 = 2d n . 



Because of the relations ( |3.4| ) and ( |3.5| ), we will not write down the results for the e 
coefficients in the Appendix A. 

We note that for the boxes (2a3) and (2a4) the Ai t *-* M-u term in ( |3.1| ) is absent as 
there does not exist a u-channel diagram for them. Furthermore, it is easy to see that 
boxes (2a3) and (2a4) go into each other with Ai t *-> M-u- For this reason we write down 
explicit results only for one of these boxes in App. A. 

Next we present the color factors and functions for the abelian type box diagram (2al). 
For this graph the sums over i in Q3.1| ) run from 1 to 7 for all the terms except for the term 
that multiplies M^, which runs from 1 to 8. We have: 

T col = - A 5 ah + (C F - ^)T b T a . (3.6) 
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ft = ln(x), / 2 = ln 2 (a;), (3.7) 



t 



h = -21n(x)ln(l + x)+21n(x)ln(-— )-2Li 2 (-x)+3C(2), / 4 = f , 



/ 5 = C(2), / 6 = ln(— f 7 = l, / 8 = ln(x)ln(l-x) + Li 2 (x). 

m 2 



dl 5 



For the nonabelian box diagram (2a2) the sums over i in (|3.1| ) run from 1 to 8 for all 
terms except for the one that multiplies , which runs from 1 to 9. For the color factor 
we obtain: 

T col = -6 ab + -Jl T b T a . (3.8) 

The relevant nine analytic functions that describe the result of evaluating the box diagram 
(2a2) are given by 

fi = / 2 = ln 2 (4)> /3 = ln 2 (4)+ 4f n- 41n (4) ln (-A) + 2 C(2), 

m 2 m 2 m 2 m z m 2 

U = hid, /b = C(2), / 6 = ln(~). /7 = 1, /s = fn, (3-9) 



h = ln(i)ln(-JU 



' m 2 m 2 ' 



In the case of the crossed box (2a4) one has nine terms for every sum in (3.1). The 
color factor for this graph takes the simple form 

Ted = V- (3.10) 

We found it convenient to define functions fi as follows: 

fi = hi, / a = fii(*-u), (3-11) 

f ?/ 

h = fdi + fdi(t^w)-f 11 -fn(t^ M ) + 21n( -)ln( 2)-4C(2), 

t V 

U = fdi, /s = C(2), /e = ln( A / 7 = ln( r 



9/5 J I \ 9 . 



f ii 
fs = ln( -)ln( -), f 9 = l. 

In order to compare with existing results we have folded the one-loop contribution 
with the LO Born term results. Using our results for the box matrix elements we have 
reproduced all the unpolarized and polarized amplitudes presented in the papers |12| and 
[|I~5H , where the Born term was folded in at the very beginning of the calculation. Note that 
after folding with the LO Born term many powers in the numerators and denominators 
cancel out, leading to a very short expressions for the box amplitudes. 

However, the above checks and the ones mentioned in the Section II did not include 
all the graphs for the gluon fusion subprocess. For a further check on the correctness of 
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the results for all of our matrix elements in Figs. 2 and 3 we have performed a global 
check on the gauge invariance by first excluding one of the heavy quark momenta from the 



independent set of momenta (see Ref. [19] for details), e.g. reexpressing p\ in terms of p\ 
and P3, then contracting the matrix elements by p±. As well, we have performed the second 
possible check, contracting our matrix elements by the momentum of the other massless 
boson p v 2 . We have verified gauge invariance for the following gauge- invariant subsets of 
diagrams: (i) When the incoming gauge bosons are photons, i.e. including graphs (2al), 
(2cl), (2c2), (2dl), (2d2), (2d3) and their u-channel counterparts, with corresponding color 
weights; (ii) For the photoproduction of heavy flavors, i.e. including all the above diagrams 
plus graphs (2a4), (2c4), (2el) and their u-channel counterparts, with corresponding color 
weights; (iii) For the hadropro duct ion of heavy flavors, which ultimately includes all the 
graphs from Figs. 2 and 3 and their relevant u-channel counterparts. We emphasize that 
the above gauge invariance checks were made separately for both color structures Cp and 
Nqi and for every existing combination of color matrices T a , T b and 5 ab , whenever they 
arise. Also, for the hadroproduction set of graphs we had to set e' = e as there is a mixing 
of UV and IR poles in the s-channel, particularly for the triangle graph (3gl) with gluons 
in the loop. 

At the end of our discussion for the gluon fusion subprocess we comment on the DRED 
result. First we note that at the one-loop level the Zip vertex renormalization constant is 
a sum of a wave function renormalization constant Z 2 , which is equal to the abelian vertex 
renormalization constant Ziy, and the term that is proportional to the Nc, which does not 
differ in DREG and DRED: 

Z 1F = Z 2 - 9 -C £ {^)N c . (3.12) 

Most importantly, the modification of the renormalization constants Z 2 and Zip in DRED 
does not affect the other renormalization constants and preserves the form of the Slavnov- 
Taylor identities for all those constants by construction. 

Next we take the set of graphs relevant to heavy flavor production in photon-photon 
collisions which are the set of graphs mentioned in item (i) above. For this abelian set 
of graphs the overall action of the renormalization counterterms amounts to completely 
removing the quark self-energy graphs with on-shell legs and adding just one term that is 
equal to one of the removed graphs, e.g. expression fl2.13| ) in DREG or ( |2.14j ) in DRED. But 



this is effectively equivalent to multiplying each of the two self-energy graphs by a factor 
of 1/2, and leaving out renormalization. Then, correspondingly, the differences between 
the two dimensional schemes for these self-energy graphs have to be halved, and simple 
counting will now show that the sum of A's for this set of graph vanishes identically. Taking 
into account the above remark on the Zip constant, e.g. separating the Cp and N c parts 
of the vertex diagrams during the renormalization procedure, one can see that the same 
statements apply to the photoproduction set of graphs given by the gauge-invariant set 
(ii) above. And the same is true also for all the self-energy and vertex graphs for the 
hadroproduction set. The remaining difference coming from the gluonic loops in graphs 
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(3gl) and (3h) cancel. This means that there is no difference between the DREG and 
DRED virtual NLO corrections for the gauge boson fusion subprocesses. 

Finally we note that the original computer output for the box diagrams was extremely 
long. The final results were cast into the above shorter form with the help of the REDUCE 



Computer Algebra System [23 



IV. ANNIHILATION OF THE QUARK- ANTIQUARK PAIR 

The graphs contributing to this subprocess are shown in Fig. 4 for the leading order 
term and in Fig. 5 for the one-loop corrections. The leading order contribution proceeds 
only through the s-channel graph. One has: 

B m = iT a T a v(p 2 )^u( Pl )u(p 3 )^v(Pi)/ s - (4- 1 ) 

Here the color matrices T a belong to different fermion lines that are connected by the gluon 
having color index a. We have again left out the factor g 2 in the above equation. In the 
Passarino-Veltman reduction for tensor integrals we use the same scalar integrals as those 
appearing in the gluon fusion subprocess, with relevant shifts and interchanges of momenta 
as needed. 

Starting again with the 2-point functions, we notice that the result for graph (5g) can 
be obtained from the one of ( [2.23] ) for graph (2h) in the gluon fusion subprocess by the 
simple replacement 

MM = M$ l) {Bir^B q9 ), (4.2) 
and all the statements after ( p.23|) are equally applicable to M( 5g ) . 



The massless quark self-energy graphs (5j) and (5k) with external legs on-shell vanish 
identically: 

M m = M m = 0. (4.3) 

The massive quark self-energy graphs (5h) and (5i) with external legs on-shell are de- 
rived analogously to the ones considered in the previous section: 

M (5h) = M (5i) = -C F B q q (i + A + ^j, (4.4) 

and the difference between the two regularizations schemes is 

A(5/i) = A(5i) = -C F B q g. (4.5) 

Results for the vertex diagrams are relatively short. Starting with graphs (5c) and (5d) 
one finds that they are proportional to the LO Born term: 

M (5c) = J B w -[-l/e / + 2/e 2 + 4/e-21n(s/m 2 )/£ + ln 2 (s/m 2 )-31n(s/m 2 )-8C(2) + 8]/6 (4.6) 
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and 

M m = 3B qq [3/e' - 4/e + \n(s/m 2 ) - 2]/2. (4.7) 
Corresponding differences between DRED and DREG results are 

A(5c) = (C F - ^)B qq , A(5d) = (4.8) 

For the other two vertex diagrams we also obtain simple expressions: 

M (5e) = {B qq [-l/e' + 3 ln(x)/3 + (1/(3 + [3)(\n(x)/e - 2Li 2 (x) - 2 ln(x) ln(l - x) 

+ ln 2 (a;)/2 - 4((2))] + 4iT a T a mv(p 2 )^ 3 u( Pl )u(p 3 )v(p 4 ) \n(x)/s 2 p}/Q (4.9) 

and 

M m = 3{B qq [3(l/e' + 2) -\n(s/m 2 )(8m 2 /s - l)/(3 2 + f M m 2 /s/? 3 ] - 4iT a T a m (4.10) 
x v{ P2 )i) Z u{ Pl )u{ P3 )v{ Pi )[\n{s/m 2 ){%m 2 /s + l)/(3 2 - 2 - 3f lld m 2 /^ 3 ]/s 2 /3 2 }/2 



with the function f lld defined in ( |2.28| ). Once again, the differences between the two regu- 
larization schemes are 

A(5e) = (C F - ^)B qq , A(5/) = ^B qq . (4.11) 

Turning to the two box diagrams we again note that because of the absence of UV poles 
there are no differences between DRED and DREG. Extensive Dirac algebra manipulations 
lead to rather compact expressions for the matrix elements. 

Since for the qq — > QQ subprocess one has two spinor "sandwiches" we cannot have 
momenta with Lorentz indices, and consequently there is no expansion of matrix elements in 
terms of Lorentz objects. We expanded the box diagrams in terms of the seven independent 
Dirac structures, the same set for each of the two box graphs. Then every Dirac structure 
is multiplied by the sums of products of a small set of analytic functions and coefficient 
functions. Thus, we have the following compact expansion for both box diagrams: 

M = iT coX {v{ P2 )Yu{ Pl )u{ P3 )^v{ PA )^2fihi (4.12) 

+ ^(P2^3w(pi)M(p3^ 1 u(p4) fi h i 

+ v( P2 )^^ a Yu{ P i)u{j)^^ a ^ u v{j>i) fihf ] 
+ mv( P2 )ft 3 u{ Pl )u{ P3 )v{ P4: ) fihi 

+ mv( P2 )Yfa^u( Pl )u( P3 )-f^ v v( Pi ) Y fib?'}- 

Note that the number of independent covariants in n ^ 4 exceeds the number of indepen- 
dent covariants in n = 4 where one has four independent covariants. 
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The sums over i in (|4.12 ) run from 1 to 4 except for the first term (it is proportional 
to the Born term), where the sum runs from 1 to 6. Below we list the color factors and 
analytic functions for the two 4-point functions of ( |4.12|) . For the graph (5a) we get: 

T coX = (T a T h )(T b T a ) ) (4.13) 

where the first parentheses in ( |4.13| ) corresponds to the summation over color indices of the 
massless fermion line and 



h = ln 2 ( — )+4f u -41n(— ) ln( j) + 2C(2), (4.14) 

t s 
f 2 = ln( -), / 3 = ln(— ), f 4 = f lld , f 5 = f u , f 6 = 1. 

m z mr 

The color factor and the corresponding functions for the second box graph (5b) are 

T col = (T a T b )(T a T b ), (4.15) 

and all the functions are obtained from the ones in ( f4.14| ) by the simple interchange t — > u, 
e.g.: 

h = ln 2 (4)+4f u (t^n)-41n(4)ln(-^) + 2C(2), (4.16) 

XL S 

f 2 = ln( -), f 3 = ln(— ), / 4 = fud, h = k(t -> u), /e = l. 

m z m z 

The coefficients hi, are given in Appendix B of this paper. However, there exists a 
partial symmetry for these box diagrams, which allows one to express most coefficients for 
the box graph (5b) through the ones of the box graph (5a). In particular, starting from 
the coefficients with superscript j > 2, we find the following general relations for all 
these coefficients: 

#[(5b)] = -h?[M](t^u), j = 2,4; i = l-M; (4.17) 
#[(5b)] = #[(5a)](t^), j = 3,5,6; 2=1-4. 

Consequently, for the graph (5b) only the coefficients hi and hf^ are presented in Ap- 
pendix B. We should also mention that all the one-loop matrix elements of this chapter 
must be multiplied by the common factor ( 2.4j) . 

At the end of this section we must again say a few words about the DRED result. In 
line of the discussion in the previous section one has a cancellation of differences between 
DREG and DRED for the heavy quark loops: the two massive self-energy diagrams (5h) 
and (5i) versus the heavy quark vertex graphs (5e) and (5f). Concerning the differences 
arising from the two massless vertex diagrams (5c), (5d) and the off-shell gluon self-energy 
(5g), these differences remain, leading to an overall difference 



A (qq _> QQ) = B qq /3. (4.18) 
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First we stress that as the above difference is proportional to the LO Born term, it is 
manifestly gauge invariant. Secondly, it can be considered as a conversion term between 
DREG and DRED for the virtual corrections for this subprocess. 



V. ABSORPTIVE PARTS 

In the previous sections we have been writing down expressions for the real parts of 
the matrix elements. However, for some physical applications as e.g. the calculation of T- 
odd observables the imaginary (or absorptive) parts of the one-loop amplitudes are needed 
as well. In addition, the NNLO calculation of heavy quark production involves also the 
product of the NLO one-loop contributions including their imaginary parts. 

In order to derive the imaginary parts of our matrix elements we have proceeded in two 
different ways. Firstly, we applied Cutkosky rules P3 to the scalar integrals that enter the 



Passarino-Veltman decomposition for our tensor integrals and calculated their imaginary 
parts directly, as discontinuities across the unitarity cut in the physical s-channel of the 
corresponding Feynman diagram. The result of these considerations may be cast into the 
form of an effective rule which states that the imaginary parts of our scalar integrals, and 
consequently the full matrix elements including real and imaginary pieces, can be obtained 
with the four simple substitutions in all of our expressions for matrix elements: 
s s 

ln( — -) — > ln( — -) — i-K, ln(x) — > ln(x) + in, (5.1) 

m 2 m 2 



ln z (—) -> hW— ) - 2z7rln(— ), \n\x) -> hf(x) + 2in\n(x). 

m m 2 m 



Note that the other logarithmic and dilog functions do not possess any imaginary parts. 
Secondly, to double-check our prescription (|5.1|), we have carefully traced the sign of the 
causal term +i0 in the propagators of the scalar integrals during the whole course of their 
derivation. We came to the conclusion that all the analytic functions are well defined in the 
physical region except for the two logarithms, which have negative arguments. They enter 
the expressions for the scalar integrals as ln(— and m(— x + iO). In other words, in all 
the expressions for our amplitudes the above two logarithms would actually appear (instead 
of ln(s/m 2 ) and ln(x)), if we had not already extracted — ir 2 terms from ln 2 (— s/m 2 — iO) 
and ln 2 (— x + iO) and considered them in the real parts of the amplitudes, and this is fully 
consistent with (|5TT|). 

In case one desires to perform crossing and obtain results for other subprocesses many of 
the analytic functions of (|3.7p , ( |3.9| ), ( |3.11|) , ( |4.14|) and ( [4.16| ) will develop imaginary parts. 



The rules for deriving such imaginary contributions were considered in detail in Ref . |25| . To 
apply those rules one needs to know the signs of the causal terms in our analytic functions. 
Thus, the only remaining thing to write down is a relative sign between the parameters of 
the analytic functions and the causal term iO for the case of our kinematics: 

s — > s + iO, t — 1 + iO, u — > u 4- iO, x — > x — iO. (5.2) 
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VI. CONCLUSIONS 



In this paper we have presented complete analytic results for the one-loop contributions 
to heavy flavor production in a closed form, including their absorptive parts. These include 
the one-loop matrix elements of the relevant partonic subprocesses ( |1 . 1| ) and ( |1.2| ) that 
are presented here for the first time 2 . We have also indicated the way of deriving matrix 
elements for the other processes that can be obtained from those presented in this paper 
by crossing to different production channels. Our results are relevant not only for various 
NLO applications, but produce that part of the next-to-next-to-leading order corrections 
to heavy flavor production corresponding to the square of the one-loop matrix elements. Of 
course, to conclude the latter task the 0(e 2 ) expansion of the relevant NLO matrix element 
is also required. We reserve this task for a future publication. 
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APPENDIX A 

Here we present the coefficients of the box contributions for the gluon fusion subprocess 
appearing in eq. Q3.1Q . 

Let us introduce the notation: 

Z\ = m 2 s — t 2 , Z2 = s + 2t, z t = 2m 2 + t, z u = 2m 2 + u, 

D = m 2 s-ut. (Al) 

First we list coefficients for the abelian type of box diagram (2al): 

h = (l/p + P)/et, b 2 = {z u -u(3)/2D, fo 3 = ZtzJtpD, b 4 = 2z u /D, 
h = (4u(3-3z u )/D, fo 6 = fo 7 = 0, fo 8 = -2(1/(3 + (3) /t; 
fog = 2z t (3/D, fog = (fog+fo 3 g)/2, bl?=st 2 (3yD 2 , 

fog = -2z t (2m 2 /t-stp 2 /D)/D, fog = -36$ - 4$, fog = 2z 2 /TD, fog = 0, 
with fog = (ztD + tz^/D 2 ; 

2 We would be happy to provide our one-loop results in REDUCE format. Please contact Z.M. by e-mail. 
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-4m 2 f3(z t + Zl /s)/D\ b$ = (fig + b^)/2, 

4m 2 t[3(D(2 + t/s) - st(3 2 )/D 3 , b$ = -Sm 2 (s(3 2 (TD + t 2 z t ) + m 2 z 2 TD/t)/tD 3 , 
-36$ - 46$, 4? = 4m 2 (2D/t + s/3 2 - 4te/s - t 2 /T)/D 2 , &g = 4z t u/stD, 
with = 4m 2 t(-^i/? 2 + uz 2 D/s 2 )/D 3 ; 

-4[3(2(2m 2 - s)z t D/s 2 [3 2 - 2m 2 s(3 2 - 3(m 2 /s + l)D + 2m 2 tz 2 /s)/D 2 , 

(&&> + 4 P 3 } )/2, 6g = -2^(t(/3 + 2m 2 ^ 2 / S 2 /?) J D + 2m 2 uz 2 [3)/D 3 , 

-A(2m 4 z 2 s(3 2 + (Am 2 TD/t 2 - z t (2m A /t - t) + 2m 2 z 2 )D)/D 3 , 

-3b$ - 46|\ = 4zt(2(D/t - z l /s + 2z 2 ) + 3tw/T - t 2 (m 2 + 2t)/T 2 )/D 2 , 

4z t (u/s -m 2 /T)/W, 

with = 2(2m 2 MZi/3 2 - (2m 2 T - 2m 2 w(w 2 + t 2 )/s 2 + t 2 (3 2 )D)/D 3 ; 

-A(3{z t D/s(3 2 + 3T 2 - t 2 (m 2 - t)/s)/D 2 , $ = (b$ + b { £)/2, 
2t{4m 2 D 2 /s 2 (3 + 2m 2 u(3D/s - t{3sT + t 2 )(3 3 )/D 3 , 
Az t ((2m 4 /t - t)D - 2stT(3 2 )/D 3 , 

-36$ - 4b ( il b$ = -8(z t /s + t 2 p 2 /D)/D, 6g = -Az t /sD, 

with b$ = 2{2t{m 2 t 2 - sT 2 ){3 2 - (Qm 2 D/s + tz t + 2m 2 t 2 z 2 /s 2 )D)/D 3 ; 

-Am 2 (3{z t + Zl /s)/D 2 , &« = (fcg + 6g>)/2, 

Am 2 tf3{{2 + t/s)/} - st/3 2 )//) 3 , &g = -8m 2 (T(D - ^),D/t 2 + s^ 2 )/^ 3 , 
-36^ - 46g, 4? = 4m 2 (2 J D/t + s/3 2 - Atu/s - t 2 /T)/D 2 , b% = 4z t u/stD, 
with = AmH(uz 2 D/s 2 - z x p 2 )/D 3 - 

-2(4m 2 /s - 3 + tz t /D)/sP, = (bff + b™)/2, 

t(4m 2 u/s 2 (3 + (2sT + t 2 )(3/D)/D, b { ™ ] = -2{2m 2 T/t - z t (2sT + t 2 )/D)/D, 
-3b^-Ab^\ b^=2(z t -2tu/s)/D, b^=0, 
with 4? ) = (2(T - m 2 w/s) - t 2 z t /D)/D; 
2(P-tz u /pD)/s, 6&) = + 6SS°)/2, 

-t 2 (4m 2 /s 2 /5 - uf3/D)/D, bff = -2(2m 2 (m 2 /t + 2) - tuz t /D)/D, 

-3bff - Abt 2 \ bff = 2{2D/s + 2s + t + t 2 /T)/D, b% } = 0, 

with 4™ } = (-2m 2 (l - t/s) + tuz t /D)/D; 

-2(2 - 2mHz 2 /sD + 3m 2 s(3 2 /D)/s(3, = + fc^)/ 2 , 

(2(T + 2m 2 t 2 /s 2 )/{3 - (2m 2 Zl - t 2 u)f3/D)/D, 

2(2m 2 (m 2 /t + 2)D - 2m 2 st[3 2 - t 2 z u )/D 2 , b { £ = -3b$ - 46^, 

-2(6m 2 + 2m 2 /s - 5mH/T)/D, bff = 0, 

with 4? = (2m 2 (l - t/s) - 2m 2 st[3 2 /D - t 2 z u /D)/D; 
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&S? = -2(2t/a - 2*Va£> + 3sTf3 2 /D)/ S (3, b { £ = (b$ + $>)/2, 

6^ = (2(-2m 2 tw/s 2 - T)//3 - (2m 2 sT - t 2 z t )(3/D)/D, 

= 2(2m 2 TD/t-2stT[3 2 -t 2 z t )/D 2 , &&> = -3b$ - Ab$ , 

= -2(3zt + 21?/s)/D, 6^ = 0, (A2) 
with 6$ = t(2m 2 /s - 1 - + tz u /D)/D; 

ci = 4/s/3, c 2 = (c 3 + c 4 /2)/2, c 3 = -(s/3 + z t //3 + 2tz 2 /sl3)/D, 

c 4 = 2(s + 3t)/D, C5 = — 4c3 — 3C4/2, C6 = C7 = 0; 

dn = 4z t /s/3L>, d 21 = (d 31 + d 4 i/2)/2, d 31 = -2T(2/sP + sP/D)/D, 

<ki = 4Tz 2 /D 2 , d 51 = -Ad 31 - 34i/2, tfci = -4/ A d 7 i = 0; 

dia = AzJsfiD, d 22 = (d 32 + d 42 /2)/2, d 32 = 2{2m 2 Zl - tuz 2 )/s[3D 2 , 

d 42 = Atz u /D 2 , d 52 = -Ad 32 - 3d 42 /2, d m = -Am 2 /TD, d 72 = 0; 

g u = 2(tp/D-4/sP), ff 2 i = (^3i+W2)/2, 3si = t{Az t D/ S p - tz 2 p)/D 2 , 

g 41 = -2(6TD-st 2 p 2 )/D 2 , g 51 = -Ag 31 - 3g 41 /2, g ei = 2tz t /TD, g 71 = 0; 

<? 12 = 4/5(1 + (m 2 ^ 2 /s - u)/s(3 2 - 2t 2 u/sD)/D, g 22 = (g 32 + g 42 /2)/2, 

g 32 = 2(-2m 2 t(QD-tz 2 )/s 2 (3 + (2m 2 s + t 2 )(3 + 2m 2 t 2 z 2 f3/D)/D 2 , 

g 42 = 4(2m 2 (m 2 - s) + z 2 - 2m 2 st 2 (3 2 /D)/D 2 , g 52 = -Ag 32 - 3g 42 /2, 

g 62 = A(Am 2 s - tu(3 + At/s) + 3m 2 t 2 /T) /D 2 , g 72 = Au/sD; 

g 13 = -A{D{m 2 + QmH/s + s)/sP + 2tUf3)/D 2 , g 23 = (g 33 + g 43 /2)/2, 

g 33 = -2((8m 2 + t)D/sP + 2m 2 t 2 z 2 /s 2 p-t(3m 2 + t 2 /s)p) + 2m 2 tuz 2 p/D)/D 2 , 

g 43 = 4(6m 4 - tz t + 2m 2 stu(3 2 / D)/ D 2 , g 53 = -Ag 33 - 3g 43 /2, 

g 63 = A(AD + 2t(2 + t/s)z 2 -3t 2 z 2 /T + t 3 z t /T 2 )/D 2 , g 73 = -4(1 - m 2 u/D)/sT; 

g u = A(m 2 D(2t/s-l)/P + 2t 3 p)/sD 2 , 92A = (g 34 + g u /2)/2, 

g 34 = 2(2(2m 2 u/s + z 2 )D/s(3-2m 2 t 2 z 2 /s 2 (3 + tz t (3 + 2t 3 z t (3/D)/D 2 , 

g 44 = 4(2T(3m 2 - s) + t 2 - 2st 2 T(3 2 /D)/D 2 , g b4 = -Ag 34 - 3g 44 /2, 

g 64 = 8(sT - 2t 2 u/s)/D 2 , g u = -At/sD; 

g 15 = A(D(m 2 z 2 /s-t)/f3-2t 2 uf3)/sD 2 , g 25 = (g 35 + g 45 /2)/2, 

g 35 = 2{2u(2m 2 - s)D/ s 2 (3 + t(2T - tu/ s)(3 + t 2 z t /s(3 + 2m 2 t 2 z 2 f3/ D)/D 2 , 

g 45 = A{6m i -2m 2 s-t 2 + 2stTu(3 2 /D)/D 2 , g 55 = -Ag 35 - 3g 45 /2, 

g 65 = A(2D -t(l + At/s)u + m 2 t 2 /T)/D 2 , g 75 =Au/sD. 

Now we list coefficients for the nonabelian box diagram (2a2): 

h = -2/et, b 2 = 0, b 3 = (s/2-t)/D, b 4 = -(z 2 /2 + z t )/(3D, 

h = -10/*, b 6 = -2z t /tT, b 7 = 2(l/e 2 + l/e + 2)/t, b s = -A/t, b 9 = A/t; 

&g) = -l/es - 2{m 2 + t 2 /s)/D, b { g = l/2s, b ( g = t 2 (z 2 /2D - 2/s)/D, 

bf = (m 2 (s-At)D/s + stz t /2-t 3 z 2 /s)/pD 2 , b$> = -4/s, 
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-8/es - 2t(2 - t/T)/D, fog = (5/e 2 + 2/e + 4)/s, fog = 8/s; 
-Am 2 (2t/D - l/sP 2 )/D, fog = 0, fo 3 g = 2z l m 2 t/D 3 1 
-2m 2 (2m 2 /s 2 f3 2 + t(m 2 (s + At) + t 2 )/D 2 )/[3D, fog = 0, 
4(m 2 (s + 3t) + mH/T)/D 2 , fog = -4u/sL>, fog = 0; 
A(2m 2 u/D + (m 2 + (3 + 4t/s)(2m 2 - s))/s(3 2 )/D, fog = 0, 
-Zt(2m 2 su/D - i)/^> 2 , 

(4m 2 (m 2 + z 2 )/s 2 (3 2 - 2m 2 / s - 2m 2 s/D + t 2 /D - 2m 4 s 2 f3 2 /D 2 )/(3D, 
0, fog = A(Dt 2 /T 2 + 2t 2 (s - t)/T + 10m 2 s + Atz t - At 2 u/s)/D 2 , 
A(2 + t/s-t/T)/D, 43=0; 

4((3m 2 + 2t)/{3 2 + 2t 3 /D)/sD, fog = 0, fog = -t(2m 2 + 3t + 2t 2 z t /D)/D 2 , 

-(3m 2 + t + (3m 2 + 2t)/{3 2 + t(m 2 s - 3t 2 )/D + 2st 4 (3 2 /D 2 )/s(3D, 

0, fog = -8t 2 z 2 /sD 2 , fog = 44/sD, fog = 0; 

-Am 2 (2t/D - l/s/3 2 )//), fog = 0, fog = 2m 2 tz 1 /D 3 , 

-2m 2 (2m 2 /s 2 (3 2 + t(m 2 z 2 + tz t )/D 2 )/(3D, fog = 0, 

4(m 2 s + 4m 2 t - t 2 + t 3 /T)/D 2 , fog = -4w/s£>, fog = 0; 

2/es + (T — t 2 /s - 4m 2 Zt/sl3 2 )/D, fog } = -1/s, 

t(z 2 /s + t(2m 2 +t/2)/D)/D, 

(2m 2 D/s{3 2 - m 2 tz 2 /s{3 2 + 2m 2 z 2 + t(2m 2 - u) + 3st 3 (3 2 /2D + t 3 z t /D)/s(3D, 



8/s, bff = 8/es + 2t(3 + At/s)/D, fog } = -3{2/e 2 + l/e + 2)/s, fog } 



2/es + (3(m 2 + iw/s) - 4m 2 te 2 /s 2 /3 2 )/ J D, b ( 22 > = -1/s, 
t(l + At/s + iu(3/2 + 2t/s)/D)/D, 

— (T - i^/s - m 2 tz 2 /s 2 (3 2 + t 2 uz 2 /2sD + m 2 st(3 2 /D)/(3D, 
8/s, fog } = 2((l/t + A/s)/e - 2(^/t - tz2/s)/D + m 2 t/TD), 
-(1/eH + 6/e 2 s - 2/et + 3/es - 2(2/4 - 3/s)), fog } = -2(1/4 + 4/s); 
-2/es + (5m 2 - 2t + t 2 /s + Az t /(3 2 + tz 2 /s(3 2 )/D, fo 2 ™ } = 1/s, 
-(2m 2 - 3t) /D + t 2 w(l/2 - 2u/s) /D 2 , 

-(Az t {sz t - 3m 2 t)/s 2 (3 2 - 6m 4 z 2 /D - m 2 st(3 2 /D - t(2m 2 /s + l))/2/3D, 
-fog) , fo& } = -fo& } - 16m 4 /T D, fog = -fo& } , 4i } = ; 
-2/es - ((8m 2 + s)z t /s(3 2 - m 2 - t 2 /s)/D, fog = 1/s, 
-(2/s - stT/D 2 - t 2 (2m 2 + t/2)/D 2 ), 
(3m 2 z t /s{3 2 - t(2T/s + 1 - (2m 2 - t/2)z t /D))/(3D, 

°51 ) °62 — °61 iDm °72 — °71 > °82 — °81 ) 
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ci = -2/sp 2 , c 2 = 0, c 3 = -s/2D, c 4 = l/2s/5 3 + z 2 /2[3D, c 5 = 0, 

c 6 = -2/T, c 7 = c 8 = 0; 

dn = -8z t /s(3 2 D, d 21 = 0, d 3 i = (-sT + t 2 )//) 2 , (A3) 

d 4 i = {z 2 /s0 2 + 2tz t /D)/(3D, d 51 = 0, d 61 = -8/D, d n = d 81 = 0; 

di2 = -8z u /s[3 2 D, d 22 = 0, d 32 = -(m 2 s + tu)/D 2 , 

d 42 = z 2 (-l/s[3 2 + 2m 2 /D)/f3D, d 52 = 0, d 62 = -8m 2 /TD, d 72 = d 82 = 0; 

#n = 2V/? 2 A 52i = 0, g 31 = -st 2 /2D 2 , 

#4i = -(2/s + m 2 z 2 /sp 2 D + stz t /2D 2 )/P, g 51 = 0, g 61 = -2t 2 /TD, g 71 = g 81 = 0; 

# i2 = -4(2 + t/s-2f//s/3 2 -3m 2 2 2 /s 2 /5 4 -2m 2 tz 2 /s/5 2 J D)/ J D, # 22 = 0, 

# 32 = s(4m 2 + ^/s + 2tV s£, )/^ 2 > 

#42 = (2(2m 4 + 2m 2 M-tT)/D + 2(4m 4 /s + m 2 + 2m 2 t/s-t)/s/5 2 + t 2 2 2 /s/? 2 L' 

- 12m 4 ^ 2 /s 3 /3 4 - 2t 3 z u /D 2 )/(3D, 

952 = 0, # 62 = 4(2D(3 + t/s) + 3te + 5m 2 t 2 /T)/ J D 2 , # 72 = 4z u /s/3 2 J D, # 82 = 0; 
013 = -4(2 + t/s + (2m 2 -u)/s{3 2 + 3m 2 z 2 /s 2 {3 4 -2tUz 2 /s[3 2 D)/D, # 23 = 0, 

#33 = s(4m 2 + tz t /s-2m 2 tu/D)/D 2 , 

#43 = (2(2m 4 - sz t -2t 2 + t 3 /s)/D + 2(4:m 4 /s + 3m 2 + I6m 2 t/s)/sp 2 -t 2 z 2 /sp 2 D 

+ 12m 4 z 2 /s 3 /3 4 - 2t 2 (tz u + z 2 u)/D 2 )/ (3D, 

953 = 0, # 63 = 4(2(2 + t/s) + 2t 2 U/TD-t 2 /T 2 )/D, # 73 = 4(t/T + z t /sf3 2 )/D, 
983 = 0; 

#14 = 4(2m 2 (s - 2u)/s 2 p 2 - 3m 2 z 2 /s 2 (3 4 + 2t 2 z t /s(3 2 D)/D, # 24 = 0, 

#34 = -(2Tu + t 2 + 2t 4 /D)/D 2 , 

g M = (3t 2 /D + (3m 2 + 18m 2 t/s -2s- 8t)/s/3 2 - 8m 4 z 1 /s 2 (3 2 D + 3m 2 z 2 /s 2 (3 4 

+ 2t 4 z 2 /sD 2 )/(3D, 

#54 = 0, g 6A = 8(2 + t/s + t 2 /D)/D, # 74 = Az t / s(3 2 D, # 84 = 0; 

#15 = 4(2m 2 z 2 /s + 3T//3 2 - Qm 2 t/s(3 2 + 2t 2 z u /D)/s(3 2 D, # 25 = 0, 

#35 = 2/£> + tz t /D 2 + 2t 3 u/ J D 3 , 

#45 = (2(2m 4 + t 2 u/s - m 2 t 2 z 2 /D) /D + t 2 z 2 /sf3 2 D - 2{2m 2 {m 2 - s) + (m 2 + s)z 2 

+ 12m 4 z t /s(3 2 )/s 2 (3 2 )/(3D, 

555 = 0, #65 = 4(-2Du/s + tz 2 -t 3 /T)/D 2 , # 75 = 4z u /s(3 2 D, # 85 = 0; 

Finally, the coefficients for the crossed box (2a4) are: 

6i = 2(2 -u/t)/s, 62 = 6i(t <->«), 63 = (2m 2 -3tu/s + s)/D, b 4 = 0, 

65 = 8(3/s - s/fu), & 6 = + ^/tT, 6 7 = 6e(t <-> u), 6 8 = -6 5 /2, 

6 9 = (l/ e 2 + l/ e + 2)s/tu; 

= -2(z u /u + 2)/s, 4? = -2(z u - 2t)/su, b$ = -2z u /su + 2t/D - tu 2 z t /sD 2 , 
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-4m 2 /tu, hi) = 4&£ ) , C = -A/eu - 2m 2 z t /TD, 

2z u (t/u + 2m 2 /U)/D, b$ = -2&£\ b$ = 2(l/e 2 + l/e + 2)/u; 

8m 2 /s 2 w, 4? = 4? = 4m 2 (2m 2 s( J D/« - t)D + tu 2 z ± ) / s 2 D 3 , 

-8m 2 /t 2 M , b<i = Ab<$, &g = -4m 2 (2(l/t - - u(3t - u)/s + m 2 t/T)/D 2 , 

-Am 2 (m 4 /U + m 2 - t 2 /s + tu/s)/D 2 , b$ = -26$, = 

4(2m 2 / M -3)/ S 2 , 6g> = 6g>, 

4? - 2(2tu/D + (4m V - tu 2 z t )/D 2 - 2m 2 tu 3 (t - u)/D 3 )/s 2 , 
-8m 2 /t 2 u, &g=4&g, 

-4(2m 2 ((l/t - 1 /s)D + u(i + 5m) / s) + tw(5m 4 + 3*u) / sT + m 2 t 2 (m 4 + tu) / sT 2 ) /D 2 , 
-4(2m 2 + 5m + 2tuz u /D)/sD, = -26$, = Am 2 z t /tTD; 

A(2m 2 t/u 2 + 4m 2 / M + 3) /s 2 , &<£ = &£\ 

&£? + 2(t(4m 2 t(t - u) - M 2 (2m 2 + 3t))D + 2m 2 t 3 u(t - u))/s 2 D 3 , 

-8m 2 /tu 2 , b$ = Ab { $, b$ = 4(2m 2 + 5t + 2tuz t /D)/sD, 

A(2m 2 ((s/u - l)D + t(u + 5*)) + ut(5m 4 + 3tu)/U + mV(m 4 + tu)/U 2 )/sD 2 , 

-2&g, 4? = -Am 2 z u /uUD; 

8m 2 (t + 2u)/s 2 u 2 , b<i = b { $ ) 

4? - Am 2 (tu(t - 2u)D + tV(t - m))/s 2 j D 3 , &g = -8m 2 /t M 2 , 
46^, 4? = 4m 2 (2m 2 s - w 2 + t 2 £//T + 2m 2 tu/T) / sD 2 , 

Am 2 (-2m 2 s(t/u + 2) + 3t 2 - M 2 (m 2 - t)/U)/sD 2 , feg = -2&£\ 4s = -IzJuD; 
-2(2m 2 + 3tu/s)/su, &£° = 4? } + 2/u, 

(-2(2m 2 s/« + t) + 2s(2m 2 w - s*)/D - t 3 u 2 /D 2 )/s 2 , 4? ) = -Am 2 /tu, 
Ab { u } + 16/u, 4? = 4/ew - 2(2D - iu)/sD, 
2/ £W + 2(2D + t 2 )/ Sj D, bff = -bff/2, bff = -3/e 2 u; 
2(-2m 2 s/u -5t-u + u 2 /t)/s 2 , b$ = 2(-2m 2 /s - 1 + 3t 2 /s 2 )/w, 
-(2s(2m 2 / M - 1) + 2t(3m 2 s + t 2 )/D - t 2 u 3 /D 2 )/s 2 , 4? = 4?\ 
4b { £\ 4? = 2{2/u - l/t)/e - 2(2m 2 (2 + u/t) + t(t - 2u)/s + t 2 /T)/D, 
2(t(-2m 2 s/u + t-2u)/s + m 2 u/U)/D, 4? = -2b { ™\ 
{I /t- 2/u)/e 2 + 2s/etu + As/tu; 

6(2t + u)/s 2 , &£? = -2(3t 2 /s 2 - l)/u, (A4) 
(2(2* + 3m) + 2t(t 2 + tu + 5u 2 )/D + ht 2 u 3 / D 2 ) / s 2 , 4? = 0, 
4bil\ b$ = -A/eu + 2(D/s + 2u 2 /s + 7m 4 /T)/D, 

2(t(2 - t/u)z u + Qm 2 u - m 2 su/U)/sD, 4? = - 2 4? , 4? = 2 (l/e 2 + l/e + 2)/u; 
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foW = g t/ / s 2 ; b (n) = _ 2{2/u _ 3t/g 2^ h {n) = ^2 + rfy p _ 3 _ ^((^ _ tu y D ^ / s ^ 

bS = 0, &£?=46&\ &£ ) = ^/et* + 2(8rn 2 -3tu/«)/A 

fog) = -2(t(2m 2 t/ M - 4m 2 + 3t) - 6m 2 M (m 2 - t)/U)/sD, b$ = -2fog\ fog = 

ci = c 2 = 0, c 3 = (u-t)/D, c 4 = c 5 = 0, c 6 = 1/T, c 7 = c 8 = c 9 = 0; 

d u = d 21 = 0, d 3 i = 2Tu/D 2 , d 41 = d 51 = 0, d 61 = 4/D, rf n = -4m 2 /UD, 

«8i = "9i = 0; 

dia = ^22 = 0, d 32 = -2tU/D 2 , d i2 = d 52 = 0, d 62 = 4m 2 /TD, d 72 = -4/D, 

d-82 — ^92 = 0; 

9u = 6/s, 5-21 = 5-11, 031 = (6 - 2t 2 /D - t 2 u 2 /D 2 )/s, ga. = 0, 951 = ^911, 

g 61 = -2m 2 t/TD, g n = -2(t + 2m 2 u/U) / D, g 81 = -2g u , g 91 = 0; 

5 12 = -12/s 2 , 522 = 512, = -2(6 - u 2 (4m 2 s + 5t 2 )/D 2 - 2t 3 u 3 /D 3 )/s 2 , 

042 = 0, 52 = 40 12 , g 62 = A(5mH/T + 3u(2m 2 -tu/s))/D 2 , 

g 72 = -4(m 2 (2t + 3u)-tu(t + 4u)/s + m 2 u 2 /U)/D 2 , g 82 = -2g 12 , g 92 = 4/D; 

013 = 012 , 023 = 012 , 033 = 032 + 4m 4 SM/D 3 , 043 = , 053 = 052, 

gm = 4(u(4D + tu)/s + 3m 2 tU/T + mH 2 /T 2 )/D 2 , g 73 = 4u(2tu - D) / sD 2 , 

083 = -20i2 , 093 = 4m 2 /TD; 

014 = 012, 024 = 012, 034 = 033 + 4t 2 UZ u /D 3 , 044 = , 054 = 052, 

6 4 = 4(3t + 4u + 2t 2 u/D)/sD, g u = 4u(-3t 2 /s + 3m 2 t/U - m 6 /U 2 )/D 2 , 

084 = -2012, 094 = 093 *~> u)] 

015 = 012 , 025 = 012 , 035 = 034 ~ 4171* st/D 3 , 45 = , 055 = 052, 

065 = 4(mH/T + w(2m 2 - 3tu/s))/D 2 , g 75 = g 65 (t <-> u), g 85 = -2g 12 , 095 = 092; 



APPENDIX B 

This Appendix contains the coefficients for the one-loop corrections to the subprocess 
qq — > QQ. For the coefficients h in ( |4.12j ) for the box diagram (5a) we have: 

h x = (t 3 /D + m 2 + 2t 2 /s)/2D, h 2 = 2(2/ es + 2/s - t/D), h 3 = -2(1 + tz t /(3 2 D)/s 1 

h 4 = -((m 2 sT + t A /s)/D - t(2T + z t j (3 2 )/ 's)/2/3D, h 5 = -4/s, h 6 = -2/e 2 s; 

= 2tT/D 2 , h ( 2 1} = 8t/sD, h ( 3 ] = -8(1 - tz t /D)/s 2 (3 2 } 

= 2z t (T/D + 2m 2 /s 2 (3 2 )/(3D; (Bl) 

= -z 1 /4D 2 , h { 2 2) = —zt/TD, hf ] = 1/D, hf = (1 - st(3 2 /D)/4(3D; 

-- t/8D, hf = 0, 4 3) = 0, hf = z t /8(3D; 

-- t 2 /D 2 , = -4/D, hf = -4z t /s(3 2 D, hf = (t Zl /D + z t /^ 2 )/s^D; 
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hf ] = st/2D 2 , hf> = 2t/TD, hf = 2z 2 /sf3 2 D ) hf = -(2z t / sf3 2 + tz 2 / D) /2/3D; 
= hf/2, hf = hf/2, hf = hf/2, hf = hT/2- 

The nontrivial coefficients for the second box diagram (5b) are: 

h x = {m 2 st/D -m 2 -4u-2u 2 /s)/2D, h 2 = -2(2(l/e + l)U/s + 1 - m 2 t/D)/U, 

h 3 = 2(1 + z t t/(3 2 D)/s, 

h 4 = (s(3 2 - 2z u - 2D/s(3 2 + Am 2 uz u /s 2 (3 2 + (m 2 z lu + sut(3 2 ) / D) /2/3D, (B2) 

hs = 4/s, h 6 = 2/e 2 s; 

hf = (m 2 s + uz u )/D 2 , hf = A(z u /U + 2u/s)/D, hf = -2(1/ (3 2 + l)(s + 2u)/sD, 

hf = ((-m 2 sf3 2 + uz lu /s)/D + (8m 4 /s + u)/s(3 2 )/(3D, 

where for convenience we have introduced the notation: 

z\ u = m 2 s — u 2 = z±(t — > u) . (B3) 

The remaining coefficients for this graph can be easily obtained from the corresponding 
ones for the graph (5a) through the relations ( f4.17| ). 
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FIGURE CAPTIONS 



Fig. 1. The t-, u- and s-shannel leading order (Born) graphs contributing to the gluon (curly 
lines) fusion amplitude. The thick solid lines correspond to the heavy quarks. 

Fig. 2. The t-channel one- loop graphs contributing to the gluon fusion amplitude. Loops 
with dotted lines represent gluon, ghost and light and heavy quarks. 

Fig. 3. The s-channel one-loop graphs contributing to the gluon fusion amplitude. Loops 
with dotted lines represent gluon, ghost and light and heavy quarks. 

Fig. 4. The lowest order Feynman diagram contributing to the subprocess qq — > QQ. The 
thick lines correspond to the heavy quarks. 

Fig. 5. The one-loop Feynman diagrams contributing to the subprocess qq — > QQ. The loop 
with dotted line represents gluon, ghost and light and heavy quarks. 
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